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ABSTRACT 

The two spectrometers on board the Infrared Space Observatory were used to observe 
the Herbig-Haro objects HH80, 81 and 80N, as well as their candidate exciting source 
IRAS 18162-2048. The fine structure lines of [Ol]63/im, [Ol]145/im and [Cll]158/im are 
detected everywhere, while [Nll]122/im and [Olll]88.3/xm are only detected toward the 
HH objects; line ratios confirm for the first time the collisionally excited HH nature 
of HH80N. No molecular line is detected in any of the observed positions. We use a 
full shock code to diagnose shock velocities Vg ~100 km s~^ toward the HH objects, 
as expected from the optical spectroscopy. Since proper motions suggest velocities in 
excess of 600 km s~^, the HH objects probably represent the interface between two fiow 
components with velocity differing by ^v^. Aside from the flow exciting source, the 
[Cll]158/im line is everywhere brighter than the [Ol]63/im line, indicating the presence of 
a Photo-Dissociation Region (PDR) all along the flow. Continuum emission from the HH 
objects and from other positions along the flow is only detected longword of ~50/im, and 
its proportionality to the [Cll]158/im line flux suggests it is PDR in origin. We propose 
that the FUV continuum irradiated by the HH objects and the jet is responsible for the 
generation of a PDR at the walls of the flow cavity. We develop a very simple model 
which strengthens the plausibility of this hypothesis. 

Subject headings: Stars: formation - (ISM:) Herbig-Haro objects - ISM: individual ob- 
jects: HH 80/81 - Infrared: ISM: lines and bands 



1. INTRODUCTION 

Stellar jets and outflows arising from low mass 
protostellar objects are considered ubiquitous, 
since essentially each protostar must go through 
an active period of mass loss to get rid of its angu- 
lar momentum to become a star (c.f. Hartmann 
19981) . These outflows are well coUimated and su- 



cited regions are detected by optical means, they 
are identifled as Herbig-Haro objects (Reipurth & 
Heathcote [19971) . 

The process of mass loss for intermediate and 
high mass protostellar objects is less well under- 
stood (Churchwell |1998| ). Observationally very 
few of these objects have collimated optical jets or 
well deflned molecular outflows (Poetzel, Mundt 
& Ray [T989t Shepherd & Churchwell pJ96|) . The 



personic and interact mostly through shocks with 
the interstellar medium. When these shock ex- 

°Based on observations with ISO, an ESA project with instruments funded by ESA Member States (especially the PI 
countries: France, Germany, the Netherlands and the United Kingdom) with the participation of ISAS and NASA. 
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Shock-Induced PDR in the HH 80/81 Flow 



Herbig-Haro 80-81 system, which is driven by a 
20 000 Lq IRAS source, is one of those rare cases 
with a well collimated supersonic jet that reaches 
flow veloc i ties o f ~ 600 - 1400 km s"^ (Marti et 
al. [19951 , 1199^ ). HH 80/81 are at the edge of 
the L291 cloud in Sagittarius (Reipurth & Gra- 
ham 119881) at an estimated kinematical distance 



of 1.7 kpc (Rodriguez et al. |1980| ). At this dis- 
tance the HH 80/81 system, with an angular size 
of ~10'.8, spans ~ 5.3 pc. The proper motion 
measurements conflrm that IRAS 18162-2048 is 
the source driving the outflow, which is part of a 
small cluster of infrared sources (Aspin & Geballe 
T992| ; Aspin et al. |199|). IRAS 18162-2048, given 



sented in §^ and discussed in 
our conclusions in 



We summarize 



its large luminosity, is likely to become a B type 
star. Most of the early work on HH 80/81 was per- 
formed at radio wavelengths, where 3.6 and 6 cm 
observations found the collimated radio jet ema- 
nating from IRAS 18162-2048, as well section 
of the counter-flow, named HH80N, which is opti- 
cally invisible (Rodrfguez& Reipurth |1989| ; Marti, 
Rodriguez & Reipurth [199^ , hereafter MR R93). 

Recently Heathcote, Reipurth & Raga ( |1998| , 
hereafter HRR) have carried a detailed analysis 
of HH 80/81 morphology and kinematics, using 
HST WFPC2 images and ground based optical 
images and spectra. They concluded in their 
study that the emission arising from HH 80/81 
is due to shocks with velocities up 600 km s~^ 
The energy release by such shocks is so high that 
thermalizes the gas to temperatures > 10^ K, i. e. 
the shocks are adiabatic. If the HH objects were 
to resemble bow shock structures (Raga & Bohm 
1986| ; Hartigan et al. |1987| ), then the optical emis- 



sion observed in Ha or [Olll] A5007 comes from 
the 'wings', while the stagnation region (the tip 
of the bow) is several times further ahead and 
invisible. This scenario, however, is derived by 
modeling the optical line proflles (HHR98); from 
the morphological viewpoint, the HH objects ap- 
pear as irregular knots that do not resemble wings 
of a bow. 

In this study we present ISO (Kessler et 
al. p.996|) spectroscopic observations of the 
HH80/81/80N system, including their exciting 
source IRAS 8162-2048. The observations and 
data analysis are presented in §H and the results 
in §0. A model to interpret the results is pre- 



2. OBSERVATIONS 

We used the two spectrometers on the ISO 
satellite to observe several locations along the 
HH 80/81 flow, including the candidate exciting 
source IRAS 18162-2048. The L ong W avelength 
Spectrometer (LWS, Clegg et al. |1996|) was used 
in its LWSOl grating mode to acquire full low res- 
olution (R~ 200) 43-197yum scans. Data were col- 
lected every 1/4 of a resolution element (equiva- 
lent to ~0.07/xm for A ^ 90/im, and to ~0.15/im 
for A ^ 90/im) with 0.5s integration time; a to- 
tal of 24 scans were accumulated, correspond- 
ing to 48s integration time per spectral element. 
LWS data processed through Off-Line Process- 
ing (OLP), version 7, have been reduced using 
the LWS Interactive Analysis^ (LIA) Version 7.2. 
The dark current and gain for each detector were 
re-estimated, and the data were recalibrated in 
wavelength, bandpass and flux. The absolute flux 
calibration for LWS in grating mode is 10% (ISO 
Handbook, iv, 4.3.2). 

The Short Wavelength Spectrometer (SWS, de 
Graauw et al. |1996| ) was used in its SWSOl "Speed 
2" grating mode. In this mode the grating is 
moved faster than the detector reset time, with 
a loss of about a factor 7 (ISO Handbook, vi 
3.3) in achievable spectral resolution, resulting 
in R~ 300. The total observing time was about 
1900s. SWS data were processed using OSIA, the 
SWS Interactive Analysis^. Dark currents and 
photometric checks were revised. The March 1998 
bandpass calibration flies have been used to pro- 
duce the flnal spectra. The accuracy of the SWS 
absolute flux calibration varies between 7% and 
35% from 2 to 40/im (ISO Handbook, vi 5.4.2). 

The flnal steps of data analysis were done using 
the ISO Spectral Analysis Package^ (ISAP) Ver- 
sion 1.6a for both LWS and SWS. Grating scans 
(LWS) and detectors spectra (SWS) were aver- 
aged using a median clipping algorithm optimised 
to flag and discard outliers mainly due to tran- 
sients. The LWS spectra were heavily fringed, and 
standard techniques available under ISAP were 
used to remove these instrumental effects. Line 
fluxes were estimated by means of gaussian flt- 



-'^LIA is available at http://www.ipac.caltech.edu/iso/lws/lia/lia.html 

^OSIA is available at http: / /www. mpe.mpg.de/www_ir/ISO /observer/ osia/osia.htm] 

^ISAP is available at http://www.ipac.caltech.edu/iso/isap/isap.htmj 
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TABLE 1 



Observations 



Object 




a(B1950) 


5(B1950) 


AOTs & TDTs 




IRAS 18162-2048 


ON 


18 16 12.9 


-20 48 49 


LWSOl-32901362 






OFFS 


18 16 12.9 


-20 50 29 


LWSOl-32901362 






OFFN 


18 16 12.9 


-20 47 09 


LWSOl-32901362 




HH80 




18 16 06.8 


-20 53 06 


LWSOl-32901350, 


SWSOl-32901351 


HH81 




18 16 07.4 


-20 52 23 


LWSOl-32901360, 


SWSOl-32901361 


HH80N 




18 16 20.7 


-20 42 53 


LWSOl-32901363 





ting. 

Table |l| lists the observed positions, which are 
also reported in Fig. |l] as the centers of the 
dashed circles (the LWS FWHM beamsize) super- 
imposed on the 6 cm VLA map of MRR93). The 
last column reports the Astronomical Observation 
Template (AOT) used, together with the 8-digits 
unique identifier of the observation (TDT). 



more intense continuum source. Secondly, molec- 
ular emission is not detected anywhere. Finally 
we note that with the exception of IRAS 18162- 
2048, the [Cll]158/im line is everywhere brighter 
than the [Ol]63/im line. The detected lines are re- 
ported in Fig. ^ for the HH objects pointings, and 
in Fig. ^ for the three central observed positions 
of the flow. 



3. RESULTS 

All objects were detected with the LWS, while 
no line or continuum is detected with the SWS 
toward HII80 and 81. Table |^ lists the measured 
fluxes for the detected lines; upper limits are given 
at the la level. The distance between HH 80 and 
81 is about one half of the LWS beamsize (~80" 
FWHM at most wavelengths, see ISO Handbook, 
IV 4.10.1) implying that each object contributes 
roughly half of its flux (lines and continuum) to 
the observed flux of its neighbour. Contamina- 
tion corrected fluxes for HH80 and 81 are listed 
just below the observed ones for each line. 

A glance at Table is sufficient to identify the 
basic results of our observations and set the guide- 
lines for the interpretation of our data. First of 
all we note that the [Nll]122/im and [Olll]88.3/im 
lines are only detected toward the HH objects; 
these observations provide the first direct ev- 
idence for the collisionally excited HH nature 
of HH80N. Although the LWS pointing toward 
IRAS 18162-2048 encompasses the radio-brightest 
portions of the thermal jet, no lines from ionised 
species (other than [Cll]158/im, but see below) are 
detected. We verified that the [Nll]122/im and 
[Olll]88.3/im lines detected toward the HH ob- 
jects would still be detectable toward IRAS 18162- 
2048, in spite of the higher photon noise from this 



3.1. The Herbig-Haro Objects 

The HH 80/81 system is the most powerful fiow 
known to emanate from an intermediate mass 
Young Stellar Object (YSO). Velocities in excess 
of 600 km s~^ are derived from proper motion 
measurements both in radio continuum (Marti, 
Rodriguez & Reipurth |1995| , |1998D and optical 
lines (HRR). In addition, line profile analysis 
of high- dispersion optical spectra suggests shock 
velocities Vg ~650 km s^^ (HRR). Interestingly 
enough however, optical spectra only reveal lines 
from species up to double ionisation state. Our 
FIR spectroscopy has the immediate advantage 
that reddening corrections, which are critical in 
the optical regime especially in this highly ex- 
tincted region (HRR), are not important. To di- 
agnose the shock conditions toward HH 80, 81 and 
80N, we compare our observations with the pre- 
dictions from plane parallel atomic shock models. 
The shock models were calculated using MAP- 
PINGS2, a code developed by Binette et al. (Bi- 
nette, Dopita, & Tuohy |1985|), w hich has been 
thoroughly tested (Pequignot |1986|) . 

A grid of models was created changing the pre- 
shock density between 10^ and 10^ cm~^, con- 
sistent with values derived by HRR from opti- 
cal lines analysis. Shock velocity were varied be- 
tween 60 and 110 km s~^; higher values, although 
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Fig. 1. — VLA 6 cm composite map of the HH 80/81 region (Marti, Rodriguez & Reipurth (1993), 
courtesy of L.F. Rodriguez); superimposed is the ISO-LWS field of view (dashed circles) at all observed 
positions. The various sources in the field are labeled. 



TABLE 2 
Observed Line Fluxes'* 



Line 

[OlJ63//m 

[Ol]145/xm 

[Olll]88.3Aim 

[Nll]122/xm 

[Cll]158//m 





IRAS 18162-2048 




HH80 


HH81 


HH80N 


ON 


OFFS 


OFFN 








104(2) 


7.1(0.5) 


8.0(0.3) 


6.4(0.2) 

3.8^ 


7.1(1.0) 

5.2'^ 


8.2(0.2) 


10.6(0.7) 


0.3(0.1) 


1.0(0.2) 


0.4(0.1) 
0.3^ 


0.32(0.04) 
0.16^ 


0.7(0.2) 


<4 


<0.3 


<0.4 


0.8(0.4) 

0.4'^ 


1.0(0.4) 

0.8^ 


1.0(0.2) 


<1 


<0.3 


<0.2 


0.35(0.08) 
O.ll'' 


0.53(0.07) 
0.47^ 


0.7(0.1) 


37(1) 


16.3(0.3) 


20.4(0.3) 


11.49(0.06) 
6.97*^ 


12.53(0.08) 
9.04^ 


17.1(0.1) 



" Units of 10^1^ W cm and la statistical uncertainties in parenthesis; upper limits are at the la level. 
Contamination-corrected fluxes. 
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Fig. 2. — Detected lines toward HH80, 81 and SON. Flux densities are normalised to 10~^^ W cm~^ 
/im-l. The dotted vertical lines represent the expected linecenter wavelength. 
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Fig. 3. — Detected lines toward IRAS 18162-2048, OFFS and OFFN. Flux densities are normalised 
to 10~^^ W cm~^ /xm-1. The dotted vertical lines represent the expected linecenter wavelength 
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Fig. 4. — [Ol]63/xm/[Ol]145/im vs [Olll]88.3/xm/[Nll]122//m diagnostic diagram obtained from the 
models of Binette et al. 1985 . Horizontal lines are of constant pre-shock density (in cm~^), while 
vertical lines are of constant shock velocity (in km s~^). 



still compatible with results from optical line ra- 
tios (HRR), fail to produce lines like [Nll]122/im 
and [Olll]88.3/im which are instead observed to- 
ward the HH objects. The preshock degree of 
ionisation and magnetic field strength were kept 
fixed at 0.1 and lO/iG respectively. These val- 
ues are not critical for the velocities we used (but 
see HRR § 8.2) and represent typical values for 
stellar jets (but see Bacciotti & Eisloffel |1999|) 
from low-mass YSOs; we will assume they are 
also valid for higher-mass systems like the one 
we are presently delaing with. This selection of 
initial values also allow us to gauge our mod- 
els with those of Hartigan, Raymond and Hart- 
mann ( 1987 ). In all instances solar abundances 
were used, since this is a reasonable assumption 
for Herbig-Haro objects (Beck-Winchatz, Bohm, 
Noriega- Crespo |1996| ). Fig. || clearly show that 
shock velocities between 90 and 110 km s~^ are 
appropriate to reproduce the emitted spectra for 
all HH objects. Pre-shock densities between 100 
and 1000 cm~^ are derived, in full agreement with 
estimates from HRR. In these conditions of low 
density J-shocks (Draine |1980| ) negligible cool- 
ing is expected from molecular lines like CO, H2 



and H2O (HoUenbach & McKee |1989|) , in agree- 
ment with our observations. Furthermore, the 
absence of FIR molecular emission also rules out 



the presence of a significant C-type shock com- 
ponent, for which H2, CO and H2O molecular 
lines would be the main coolants (e.g., Kaufman 
& NeufeldjOTp. 

These low velocity shocks plus the width of 
the observed optical lines suggested to HHR that 
the HH 80/81 objects correspond to "wings" of a 
highly adiabatic bow shock that strikes the sur- 
rounding gas at ~ 650km s~^; in such wings the 
shock front would be oblique with respect to the 
direction of motion, resulting in lower shock ve- 
locities. A shock velocity of ~ 100 km s~^ and 
flow velocities (from the proper motions) as high 
as 600 km s~^ can be also reconciled if one consid- 
ers that the mass loss from the protostar is time 
dependent, and that the observed configuration 
of HH 80/81 was preceded by previous ejection 
events. If this is the case, the new ejected gas finds 
the circumstellar gas already in motion, and it is 
the interaction between these two flows that can 
create relative shock velocities of ~ 100 km s~^ 
so (Raga et al. 



or 



1990 



Stone & Norman, |1993| ). 
Multi-epoch VLA observations (Marti, Rodriguez 
& Reipurth [T995| , [1998D of the jet's radio knots 



provide supporting evidence for episodic mass loss 
from IRAS 18162-2048. 

A similar process has been invoked to under- 
stand the flow characteristics of the HH 1/2 sys- 
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tern, where the proper motions are as high as 
450km s~^ in both atomic and molecular H2 gas 
(Herbig & Jones |1981| ; Rodriguez et al. |199CI| ; 
Eisloffel, Mundt & Bohm 1994 ; Noriega-Crespo 
et al. |1997| ), but the shock themselves are ~ 
150 - 200 km s'^ (Har tigan et al. [1^571 ; Noriega- 
Crespo, Bohm & Raga pJ89D . In the HH 1/2 out- 
flow, however, there is clear evidence of a previous 
outburst event marked by the presence of an older 
bow shock structure ~ 10 times farther away from 



In a PDR, electrons are released in the gas phase 
by FUV-irradiated dust grains via the photoelec- 
tric effect and heat the gas (Tielens & HoUen- 
bach 1985). Theory and observations generally 



central outflow source (Ogura p,995|) 



3.2. The Jet 

The LWS pointing toward IRAS 18162-2048 en- 
compasses the base of the radio thermal jet, which 
is the strongest radio emitter in the flow. Yet no 
[Olll]88.3/im or [Nll]122/im emission is detected, 
indicating that the ionised material of the jet is 
at a lower temperature compared to IIII80, 81 
and 80N, or that the ionisation does not result 
from shocks. Interestingly, the radio spectral in- 
dex of the jet is consistent with that of an ionised 
wind, but signiflcantly differs from that of the 
HH objects which instead manifest a possible syn- 
chrotron component (MRR93). This behavior is 
different from that of the HH 1/2 system, for in- 
stance, where the central source VLA 1 (Pravdo 



et al. |1985| ) has a similar positive radio spectral 
index to that of HH 80/81 jet, but for the HH ob- 
jects themselves it is flat, indicating optically thin 
free-free emission (Rodriguez et al. 1990| ). 



3.3. The PDR 

In the shock conditions diagnosed by our FIR 
spectroscopy (see Sect. |0| ), [Cll]158/im con- 
tributes about 20% of the [Ol]63/xm cooling. Since 
[Cll]158/im is everywhere brighter than [Ol]63/im 
(except toward IRAS 18162-2048) the [Cll]158/im 
emission must have another origin, namely PDR. 
Evidence for a PDR-like emission is also of- 
fered by the FIR continua detected toward the 
HH object and the OFF positions (the latter 
cannot be entirely justifled with contamination 
by the stronger continuum source IRAS 18162- 
2048). Fig. ^ shows the SEDs observed toward all 
pointed positions with the LWS. We note that the 
continuum observed toward the two OFF posi- 
tions cannot be justifled with contamination from 
IRAS 18162-2048. The continua from HH80 and 
81 have been corrected for reciprocal contamina- 
tion (see enf of §0), and detector spectra have 
been stitched for cosmetic purposes. 



agree on the fact that between 0.1% and 1% of 
the incident FUV flux is converted into gas heat- 
ing via this mechanism, and subsequently released 
mainly via "cooling" of the [Cll]158/im line. The 
rest of the incident fleld is absorbed by dust and 
reprocessed to FIR wavelengths. We see from 
Fig. 1^ that the observed SEDs contain the bulk 
of FIR emission, so that their integral is a good 
representation of the reprocessed FUV fleld; the 
ratio of observed [Cll]158yum to the integrated 
SEDs yields values 0.002< Xc <0.013, in very 
good agreement with expectations for a PDR. The 
only obvious exception is IRAS 18162-2048, whose 
SED is dominated by radiation from the central 
YSO. 

Assuming that the pre-shock densities diag- 
nosed for the HH objects (Sect. |3?T|) represent 
the average conditions surrounding the HH 80/81 
flow, we can use the observed [Cll]158/im flux 
to estimate the intensity of the irradiating FUV 
fleld. Using the Web Infrared Tool Shed (Wolflre 
et al.) and assuming complete flUing of the LWS 
beam with fllling factor of unity we estimate Go 
-200, 400, 30 000, 200, 40 and 30 for the 6 LWS 
pointings from north to south (see Fig. |l]), where 
Go is the FUV fleld intensity expressed in units of 
1.6x10"^ ergs cm~^ s~^ (Habing |1968| ). In these 
conditions the PDR emission could account for 
most of the [Ol]63/im and [Ol]145/im fluxes every- 
where along the flow, including the HH objects. 
In this case the [Ol]63/im/[Ol]145/im ratios would 
imply densities n < 10^ cm~^. 

We searched several degrees around the 
HH 80/81 area for OB stars which could be re- 
sponsible for this irradiation level, flnding none. 
An obvious candidate is of course IRAS 18 162- 
2048 itself. Assuming a BO ZAMS spectral class 
(MRR93), the resulting FUV fleld at the vari- 
ous LWS pointings would be Go ~450 at OFF- 
N and OFF-S, and ~60 at the HH objects; the 
agreement with the Go values estimated from the 
[Cll]158/im line can be considered satisfactory. K- 
band images of IRAS 18162-2048, however, clearly 
resolve the IRAS source in a small cluster of at 
least 3 sources (Aspin et al. |1994| ); redistribution 
of the total bolometric luminosity of the IRAS 
source among the cluster members, adopting the 
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Fig. 5. — Spectral energy distributions observed with LWS toward all pointings. The SED from IRAS 
18162-2048 has been divided by 5 to fit it into the plot range. The SEDs of IIII80 and 81 have been 
corrected for reciprocal contamination. 



IMF from Miller & Scalo (|T97|), would assign 
only half of this luminosity to the most massive 
member of the cloud (e.g., Molinari et al. |2000| , 
and references therein). This decreases the Ly- 
man continuum photon flux by more than one or- 



der of magnitude and reduces Go by nearly two 
orders of magnitude, ruling out IRAS 18162-2048 
as the PDR illuminating source. 

An intriguing possibility is that the illuminat- 
ing sources for the observed PDR emission are 
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Fig. 6. — Plot of the [Cll]158/im line flux vs the 6 cm continuum flux from the knots of the jet 
encompassed by the LWS beam. When errorbars are not visible is because they are smaller than the 
symbols. 
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the HH objects and the jet. In Fig. ^ we report 
the [Cll]158/im hne flux as a function of the 6 cm 
flux measured by MRR93. The filled symbols rep- 
resent the HH objects and IRAS 18162-2048; for 
each object we consider the radio knots encom- 
passed by the LWS beam at each position. The 
correlation between Fc, and Fgcm would seem to 
suggest that the ionised material and the amount 
of cooling by the PDR maybe somehow related. 
We note that the two OFF positions do not fit 
with the correlation in Fig. In the next para- 
graph we will propose a semi-empirical model to 
verify and quantify the Radio-PDR connection. 



4. A MODEL FOR THE SHOCK/JET-PDR 
CONNECTION 

The influence of shock-originated radiation on 
the surrounding medium in HH flows is not clear. 
NH3 and HCO"'" enhancements have been de- 
tected in the vicinity of HH80N, HHl and HH2 
(Girartet al. |T^; D avis, D e nt fc B ell Bur- 
[I990t Torrelles et al. [T9U2| , [T993D and it 



nell [19901 Torrelles et al. 
has been suggested that the UV fleld generated 
in the HH shocks may be responsible. Wolflre 
& Konigl (|1993| ), and more recently Raga & 
Williams ( pOOOj ), have shown that the chemistry 
of blobs in the vicinity of HH objects can be influ- 
enced by the passage of the HH object, although 
the predicted morphologies do not match the ob- 
servations. 

We propose a scenario where the ionised mate- 
rial which is recombining behind the shock front 
and in the jet, and which emits free-free radio con- 
tinuum, is also responsible for a UV fleld which 
illuminates internally the walls of the flow cavity; 
a PDR is there produced, which mainly cools via 
the [Cll]158/im line. If this is correct, then we 
should be able to express the two observables, i.e. 
the radio continuum flux and the [Cll]158/im line 
flux, function of a set of physical parameters 
which characterize the HH shocks, the jet and the 
PDR. 

Radio emission from shocked material has been 



modeled by Curiel, Canto and Rodriguez (|1987|) ; 
in the optically thin regime, the free-free emission 
from the recombination region can be written as 



^shv 



1.8410"^ r 



n -0.1 



10 GHz 



T, 



0.45 



n. 



010 



where 6 is the angular diameter in seconds of 
arc of the recombination region, T4 is the elec- 
trons temperature in units of lO^K, no^^ is the 
pre-shock density in units of 10 cm~^, and Vg^ is 
the shock velocity in units of 100 km s~^. This 
model applies to HH80, 80N and 81, whose high 
ionization lines clearly trace a high velocity shock. 
In the absence of similar evidence toward the cen- 
tral source and the two OFF positions (Tab. |), 
where the LWS beam encompasses most of the jet 
radio knots, we will assume that the ionisation 
does not results from shocks. Reynolds (|1986|) 



developed a model to predict the properties of ra- 
dio continuum emission from a collimated, ionised 
thermal jet. His expression of the radio flux at 
any speciflc frequency is a complicated function 
of jet parameters like the collimation, the density, 
temperature and ionization radial gradients, the 
initial distance where the jet is injected, the jet's 
width at its base. Because we cannot indepen- 
dently flx any of these parameters, the diagnostic 
power of this model will be very limited in the 
present case. 

Let us now quantify the energy released by the 
shock recombination region, in the portion of the 
UV continuum, between 912A and 2066 A(or be- 
tween 13.6eV and 6eV), which is effective in PDR 
illumination. It is widely accepted (Dopita, Bi- 
nette & Schwartz |1982| ) that the dominant contri- 
bution to UV continuum in this wavelength range 
comes from hydrogenic 2s— s-ls two-photon decay. 
Theoretical models (Shull & McKee 1979| , also 
conflrmed by our grid of computed shock mod- 
els (§p.l[)) suggest that two-photon emission from 
He° and He+ can be neglected in comparison. The 
H° two-photon emissivity is given by 



1 /iz/ _i _-. 

3p = Pp/voN2s erg s sr cm 

47r uq 



Hz- 



where is the Lya frequency, ^u/uq is the 
probability, symmetrical around ^vq, that a pho- 
ton is emitted with frequency f/vo and for which 
we adopted the analytical approximation of Nuss- 
baumer & Schmutz ( 1984 ). is the population 
density of level 2s and can be found equating the 
recombination rate from higher states to 2s, with 
the total 2s— *ls two-photon decay rate. We ob- 
tain (Emerson |1996| ): 



Vs,[l + 3.483W,, - 2.745] mJy 



N2s ~ 0.3 a'.^^NeNiA:,^ cm'^ 



(3) 
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Fig. 7. — Sketch of the proposed scenario for the HH-PDR connection. The HH object is the grey 
knot at the origin of the coordinate system, and the cyhnder represents the cavity excavated by the 
flow. The portion of big circle at the top of the figure indicates the ISO-LWS beam. 



where Ne and Nj are the electron and ion num- 
ber densities, and ~ 8.23 s~^ is the two- 
photon transition probability. The total recom- 
bination coefficient to all hydrogen excited levels 
can be written as (Hummer & Seaton |1963| ): 



1.62710"^^ T4 [1 - 1.657 1.0^10^4 + 

(4) 



0.584 T^/^ ] cm^ s"^ 



where T4 is T/10^ K. Using Eqs. (|), (|) and 
(^, we are able to write the power radiated by 
the recombination region in the FUV as: 



Lfuv = 47r Vrec Vce / judu crg s"^ (5) 

JGeV 

Vrec is the volume of the recombination region, 
while r]ce is a factor which accounts for a collisional 
enhancement of the 2s level population above the 
values predicted by pure recombination. Such an 



enhancement can be determined from a compar- 
ison of the predicted two-photon spectrum with 
the observed UV continuum from HH objects. 
Dopita, Binette & Schwartz ( |1982| ) find values 
2.8 < rjce < 13.3; interestingly, rjce is found to 
be inversely proportional to the degree of excita- 
tion of the HH object as measured, e.g., from the 
[Olll]5007A/[Ol]6300A ratio. Eq. (§ makes the 
implicit assumption that the ionised clump is op- 
tically thin to the FUV radiation. It is easy to 
show that for typical dust grains of radius O.l/xm 
and density of 3 gcm~^, the optical depth along 
the diameter of a spherical clump of radius r and 
particle density n can be written as: 



r,. 



l.SlO^^Q^nrd (6) 

where n is in cm~^, r is in seconds of arc and 
the distance d is in parsecs. In the FUV range 
the absorption coefficient ~ 1 (Draine & Lee 
1984] ), and for the typical parameters that we will 
derive (see Table 0) we obtain tfuv ~ 0.8. We 
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emphasize that this number refers to the longest 
path across the clump, so only a small portion at 
the far side of the clump, with respect to any line 
of sight, will be only partially thick. 

The FUV field emitted by the HH object is in- 
tercepted by the flow cavity walls, and a PDR is 
there formed. A fraction Xc of the incident FUV 
field is predicted by PDR models (Tielens & Hol- 
lenbach 1985| ) to be reradiated via the [Cll]158/im 
line; our observations allow us to determine (§ |3.3D 
values of Xc in very good agreement with model 
predictions. Since we want to compare the pre- 
dicted [Cll]158/im line flux radiated by this PDR 
with our ISO-LWS observations, we are interested 
in the portion of the flow cavity which is encom- 
passed by the instrument field of view. In other 
words, we need to estimate the fraction fc of the 
^Fuv emitted by the HH object which is inter- 
cepted by the flow cavity in our LWS beam. The 
situation is sketched in Fig. |^, which serves as ref- 
erence for the following discussion. R is the radius 
of the cavity, and Rf, is the radius of the ISO-LWS 
beam. For simplicity, we will assume the HH ob- 
ject lying at the center of our coordinate system 
as pointlike with respect to R. The element solid 
angle at the cavity wall as seen from the origin is 



R dO dz cos 



(7) 



where the cos0 factor accounts for the projec- 
tion, perpendicular to the line of sight from the 
origin, of the element area of the cavity wall. Ex- 
pressing z and r as functions of R and 0, we obtain 



dn 



cos 



de 



The total solid angle VL^ under which the HH ob- 
ject sees the internal cavity walls (from z=-Rb to 
z=Rf,) is obtained by integrating Eq. (H) over 27r 
in d6', and from -arctan(Rfe/R) to arctan(Rf,/R) 
in d0. The FUV field is emitted isotropically, so 
that 



Jc = -r = sm 
47r 



, Rb 

arctan — 
R 



(9) 



We can finally write the predicted flux of the 
[Cll]158/im line from a system at the distance D 
as: 



10 Lpuv fc Xc 



Wcm 



(10) 



5. DISCUSSION 
5.1. The Herbig-Haro Objects 

To understand which are the critical parame- 
ters in our model, it is useful to show in Fig. | 
a diagnostic diagram that presents the relation- 
ship between our two observable quantities, Fgcm 
and Fcii, for various parameters sets. We will first 
consider the case of the HH objects, where Eq. (|l]) 
holds for the emitted radio continuum flux. For 
each grid in Fig. ^, the horizontal lines are sites 
of constant density (n in cm~^) while the verti- 
cal lines indicates constant TZ=n/no compression 
ratio (density over pre-shock density). The full- 
line grid is computed for Vs=100 km s~^, rice='2, 
Xc=0.007 and thh=^" ■ Decreasing the shock ve- 
locity to 70 km s^^shifts the grid to lower radio 
fluxes (dashed), while a lower FUV— >-Cn, Xc, con- 
version factor brings it down in [Cll]158/im flux 
(dash-dotted, the opposite effect is obtained by 
raising the r/ce enhancement factor). Decreasing 
the radius of the HH object influences both quan- 
tities (dotted). We find no appreciable effect from 
the particular choice of the cavity radius R. It 
is reassuring that for reasonable choices of the 
model parameters we can reproduce the 6cm and 
[Cll]158/im fluxes observed toward the HH ob- 
jects. The number of free parameters in our model 
can be significantly decreased thanks to the avail- 
able observational evidence. The value of Xc for 
each object has been determined using our FIR 
continuum and [Cll]158/im measurements (§ j3.3| ). 
The radius of the HH objects is taken from the 
radio maps of MRR93; the objects are clearly re- 
solved, and the full size (after a tentative beam 
deconvolution) is ~10". The value of rjce is com- 
puted using its relationship (Dopita, Binette & 
Schwartz |1982|) with the degree of excitation mea- 
sured by the [Olll]5007A/[Ol]6300A ratio; we use 
the [Olll]5007A and [Ol]6300A line fluxes from 
HRR to derive values of 2 for the HH objects 
and 13 for the jet. We adopt shock velocities of 
100 km s~^ for all HH objects from the measured 
[Olll]88.3/im/[Nll]122/xm ratio. Model results are 
weakly dependent on the cavity radius, for which 
we adopt a fiducial value of R=10". We can then 
vary the plasma density and TZ to fit the data, 
and the results are in Cols. 8-10 of Table |^. 

For all HH objects we obtain densities which are 
below lO^cm"^, for which 2s 2p collisional pop- 
ulation would depress the two-photon continuum 
due to Lya decay. The presence of a transverse 
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TABLE 3 



Model inputs and results for HH objects 



Object 


Xc 


thh 


R 


Vs 




FUV 


n 


7^ 






(") 


(") 


(km s"""^) 




(Go) 


(cm~'^) 




HH80 


0.011 


5 


10 


100 


2 


170 


4200 


13 


HH81 


0.013 


5 


10 


100 


2 


180 


4500 


10 


HH80N 


0.007 


5 


10 


100 


2 


670 


8400 


13 



(Bo 1-Vs) magnetic field limits the compression 
ratio; we can write (HoUenbach & McKee |1979| ) 



5, 



76.7 



100km s 



-1 



:iii 



Using the parameters from Table ^ we see that 
i?o± varies from 100 to 200 /xG. We ran additional 
shock models using these large Bo± values and 
found that the only way to reproduce the observed 
line ratios is to assume a pre-shock ionization frac- 
tion > 0.5. This is higher than the higher values 
(0.2-0.3) found toward several HH-jets by Bac- 



ciotti & Eisloffel ( p.999|) ; on the other hand the 
HH 80/81 fiow is the most powerful in its class 
and may be peculiar in its ionization fraction as 
well. 

5.2. The Jet 

As far as the other positions along the jet are 
concerned, our model has less diagnostic power 
because of the high number of free parameters 
in the Reynolds ( |1986|) model for radio emission 
from ionised jets. 

Fig. ^ is the analogous for the jet of Fig. ^, 
where grids of models are plotted function 
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0.01 0.10 1.00 10.00 

Fig. 8. — [Cll]158/im line fiux as a function of 6 cm radio continuum fiux (from Eq. |ip for different 
choices of HH model parameters. Each grid shows the locus of points of constant density n (horizontal 
lines) and constant 7^=n/no compression ratio. The full-line grid is computed for fs=100 km s~^, 77ce=2, 
Xc = 0.007 and thh=^" ■ The other grids have fs=70km s~^ (dashed), Xc = 0.002 (dash-dotted) and 
yhh='^-5" (dotted). The full circles represent the observations. 
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Fig. 9. — [Cll]158/im line flux as a function of 6 cm radio continuum flux (from Reynolds 1986) for 
different choices of jet model parameters. Each grid shows the locus of points of constant density n 
and jet's initial half-width in arcseconds. The full-line grid is computed for [e, q„, q^, qa;] = [0.45, -0.9, 
-0.6, -0.5]. The other grids have qr=-l (dashed), e = 0.2 (dash-dotted) and q„=-2 (dotted). The circles 
represent the observations; the open symbols mean an uncertain assignment of the radio flux. 



of density n in cm~^ and half-width wq, in arcsec- 
onds, of the jets at its base. While Xc and rj^e are 
kept fixed to 0.002 (from our observations) and 
13 (see above), the jet parameters are essentially 
free. The full-line grid represent the "pressure- 
confined" jet in Reynolds terminology, where the 
density, temperature and ionization radial gradi- 
ents of the jet have values [q„, qr, qa;] = [-0.9, - 
0.6, -0.5]; in this particular jet the degree of col- 
limation, expressed by a number e which controls 
the jet's width function of radial distance 
(woe r*^), is e = 0.45. Another critical param- 
eter, which we cannot independently fix, is the 
distance from the central source where the jet is 
injected; the full-line grid in Fig. ^ is obtained for 
ro=10 -R0. This model predicts too much radio 
flux compared to the observations; other classes 
of standard models in Reynolds study are even 
less successful. Although of uncertain physical 
meaning, other combinations of the jet parame- 
ters can bring the predicted radio fluxes closer 
to the values observed along the HH 80/81 jet. 
Steepening the density gradient to qn=-2 (dotted 
grid) will decrease the radio flux by more than 
an order of magnitude; a less pronounced effect 



in the same direction is obtained by decreasing 
the temperature gradient to, e.g., qT=-l (dashed 
grid). Providing additional coUimation by low- 
ering e down to, e.g., 0.2 (dashed-dotted grid), 
will also contribute to a better match between 
the model and the data. Finally, decreasing the 
jet-injection distance from the source, ro, will also 
work toward lowering the predicted radio flux. Al- 
though there is no single set of jet parameters 
that can fit the data, our model suggests steep 
density and/or temperature radial gradients along 
the jet, as well as a very high degree of coUima- 
tion which is not ruled out by the observations 
(MRR93). Looking at the shape of the models 
grids in Fig. ^, it seems that a single jet model, 
i.e. a single model grid, cannot reproduce the 
[Fcii, Fgcm] values of IRAS 18162-2048 and the 
two OFF positions. This is not due to the mis- 
alignment between the observed OFF positions 
and the jet axis. A better alignment would not 
have changed the radio flux (which we arbitrar- 
ily assigned based on the 6 cm map of MRR93); 
likewise, a half-LWS beam shift in the ISO ob- 
served positions would certainly not result in a 
one order of magnitude decrease of the [Cll]158/im 
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hne fluxes, as required to bring the OFF positions 
onto the model grid in Fig. ||. Rather, it is plau- 
sible that the central (brighter) portion of the jet 
considerably contributes to the FUV irradiation 
of the flow cavities at the two OFF positions. 
This would indeed correspond to increasing the 
assigned radio flux to the two OFF positions in 
Fig. |. 

5.3. The Fcn-FQcm relationship 

Although it was one of the motivations to de- 
velop our model, one of the consequences of the 
model is that a Fcn-Fgcm relationship is difficult 
to justify. The radio flux has a different origin 
in the HH objects and the jet, there is a priori 
no reason why IRAS 18162-2048 (where we model 
the radio flux as coming from a jet) should line up 
with the three HH objects (where we model the 
radio flux as coming from the post-shock region) 
in Fig. |. The relative position of the HH objects 
themselves is also a priori depending on a high 
number of free parameters. Based on four points 
only, we must then conclude that the Fcu-Fecm 
relationship in Fig. |^ appears to be fortuitous. 
However, we will check for similar occurrences in 
other HH/jet systems where [Cll]158/xm data are 
available and where the densities are low enough 
[n < 10^ cm~^) that a shock/jet <-i>PDR connec- 
tion can be expected. 

6. CONCLUSIONS 

We have performed a far-IR spectroscopic study 
of the HH 80/81 system. Line ratio analysis con- 
firms for the first time the Herbig-Haro nature 
of the nebulosity HH 80N, which probably repre- 
sents the head of the counterflow to HH 80 and 81. 
We reveal shock velocities of the order of 100 km 



s~^ in correspondence with the HH objects, while 
lower excitation conditions appear to be present 
elsewhere along the radio jet. A comparison with 
proper motion velocities in excess of 600 km s~^ 
indicate that the shocks arise at the interface be- 
tween two fast-moving flows. 

Besides shock-excited emission, an important 
PDR contribution is present all along the bipo- 
lar flow, where densities below 10^ cm~^ are also 
diagnosed. Using a simple model, we have pro- 
vided quantitative arguments supporting the idea 
that the FUV field radiated by the ionised mate- 
rial of the recombination regions in the HH ob- 
jects and of the jet emanating from IRAS 18162- 
2048, is able to induce the formation of a PDR 
in the immediately surrounding medium (i.e. the 
flow cavity walls). This would provide further ev- 
idence that the jet/HH own radiation field affects 
its surrounding medium in a measurable way. At- 
tempts to model the chemistry of outflows should 
not ignore the influence of the radiation field of 
the shocks responsible for the acceleration of the 
outflow itself. 
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